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Muon-spin rotation (/iSR) experiments are often used to study the magnetic field distribution in 
type-II superconductors in the vortex state. Based on the determination of the magnetic penetration 
depth it is frequently speculated — also controversially — about the order-parameter symmetry of the 
studied superconductors. This article reports on a combined /iSR and magnetization study of the 
mixed state in the cuprate high-temperature superconductor Lai.83Sro.i7Cu04 in a low magnetic 
field of 20 mT applied along the c axis of a single crystal. The macroscopic magnetization measure- 
ments reveal substantial differences for various cooling procedures. Yet, indicated changes in the 
vortex dynamics between different temperature regions as well as the results of the microscopic /iSR 
experiments are virtually independent of the employed cooling cycles. Additionally, it is found that 
the mean magnetic flux density, locally probed by the muons, strongly increases at low tempera- 
tures. This can possibly be explained by a non-random sampling of the spatial field distribution of 
the vortex lattice in this cuprate superconductor caused by intensified vortex pinning. 

PACS numbers: 74.25.Ha, 74.25.Wx, 74.72. Gh, 76.75.+i 



I. INTRODUCTION 

The muon-spin rotation (^SR) technique is a power- 
ful tool to study the local magnetic field distribution 
in solids It is also successfully employed to gain valu- 
able information on type-II superconductors by probing 
the magnetic field distribution P{B) generated by a vor- 
tex lattice (for a revievif mostly regarding cuprate high- 
temperature superconductors, see, e. g. Ref. i). If the 
vortex-lattice configuration is ordered and static, from 
the obtained P{B) the characteristic length scales of the 
superconductor like the magnetic penetration depth A or 
with reservations the Ginzburg-Landau coherence length 
^ can be extracted in a reliable way. However, the vor- 
tex arrangement in cuprate superconductors is commonly 
not an ideal static two-dimensional hexagonal flux-line 
lattice. Individual vortices are always subject to dis- 
placements due to pinning, e. g. at grain boundaries or 
lattice defects. Also either thermal or quantum fluctu- 
ations introduce dynamics and may cause a reordering 
and relaxation of a non-equilibrium vortex lattice. Re- 
views on this subject can be found, e. g. in Refs. iandi. 
The presence of those imperfections changes the ^SR field 
distributions: while weak random pinning leads to a sym- 
metric broadening of P{B),^ thermal fiuctuations might 
even change the observed asymmetry of the lineshape as 
seen in the example of the vortex-lattice melting in the 
Bi2Sr2CaCu208+5 compound.^ Moreover, the lineshape 
asymmetry has also been found to depend on three-body 
correlations between the vortices as seen in the so-called 
"vortex-glass state" in Lai.9Sro.iCu04.^ 
Therefore, in order to obtain accurate information on 
the superconducting-state parameters, it is preferable 
to complement transverse-field (TF) fiSK measurements 



in the mixed state by investigations using other exper- 
imental techniques. Here, we report on combined TF 
/iSR and magnetization studies of the mixed state of a 
Lai.83Sro.i7Cu04 single crystal in a low magnetic field of 
/io-ff = 20 mT applied parallel to the c axis of the crys- 
tal. A detailed analysis shows that although the mea- 
surements are not susceptible to small variations in the 
order-parameter symmetry, the data can be overall con- 
sistently described by taking into account a single energy 
gap in the quasiparticle excitation spectrum with dx2-y2 
symmetry and a change in the sampling of the spatial 
field distributions by the muons due to vortex-pinning 
effects. 



II. EXPERIMENTAL DETAILS 

A cylindrical Lai.83Sro.i7Cu04 single crystal with a di- 
ameter of 5.5 mm and a height of 10 mm was used for the 
fiSK studies. The crystal was grown by the travelling- 
solvent fioating-zone technique similar to the one re- 
ported in Ref. The subsequent characterization by 
Laue X-ray diffraction showed the c axis pointing perpen- 
dicular to the cylinder axis. X-ray powder diffraction on 
a sample of the same batch indicated a single crystalline 
phase. The magnetization studies have been carried out 
in a Quantum Design 5T MPMS SQUID magnetometer 
on a small piece of approximate dimensions 5x4x2 mm^ 
of the same batch with the c axis along the shortest edge 
of the sample. The critical temperature Tc is about 37 K 
(cf. Fig.[T|). The magnetization measurements in a mag- 
netic field in the range from a few millitesla up to several 
tesla applied parallel to the c axis of the crystal were 
performed using different field-cooling (FC) procedures: 
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(i) fast cooling from above Tc to T = 5 K with a high 
rate of w — 20 K/min and measuring during warming up 
the sample (FFCW), (ii) slow cooling from above Tc to 
T = 5 K with a low rate of « —0.7 K/min and measuring 
during this cooling (SFCC), and (iii) measuring during 
warming up from the slowly field-cooled state obtained 
through (ii) (SFCW). 

The /xSR experiments with magnetic fields up to 20 mT 
applied along the c axis were performed at the GPS spec- 
trometer at the 7rM3 beam line at the Swiss Muon Source 
at the Paul Scherrer Institute)^ The used helium flow 
cryostat allows to cool the sample from T = 50 K to 
T = 1.6 K in about two minutes at the maximum cooling 
rate. It is therefore possible to do the /zSR measurements 
using equivalent cooling procedures as for the magneti- 
zation studies. It should also be noted that the ratio of 
the sample dimensions parallel and perpendicular to the 
magnetic field has been similar for the magnetization and 
the /iSR measurements. Thus, the demagnetization ef- 
fects on the respective samples are similar for both stud- 
ies and overall small due to the small absolute value of 
the susceptibility in the field-cooled mixed state. 
For the /iSR experiments spin-polarized positively 
charged muons with an energy of ^ 4 MeV are implanted 
into the sample and thermalize. The spin ensemble in- 
teracts with its local environment until the muons decay 
(r^ = 2.197 lis) and emit the formed positrons prefer- 
entially in the direction of the muon spin at the time 
of decay. Thus, by detecting the muons at their im- 
plantation time and the positrons after the decay the 
temporal evolution of the muon-spin polarization in a 
sample may be recorded and by that information about 
the local environment of the muons is obtained. In a 
static local magnetic field Bioc with a non-zero compo- 
nent perpendicular to the muon spins, these undergo a 
Larmor precession with a frequency uj = ^^B\oc, where 
7^ = 27r X 135.54 MHz/T is the gyromagnetic ratio of 
the muon. In a homogeneous transverse field the mea- 
sured field distribution would therefore ideally be a Dirac 
5 function, whereas P{B) has a finite width when spa- 
tially inhomogeneous fields are probed by the muons. In 
case the local fields within the mixed state of a type-II 
superconductor are sampled, this leads to a A- and ^- 
dependent characteristic asymmetric P{B) with a "high- 
field tail" originating from contributions of muons stop- 
ping in and close to the vortex cores.^^ Further details 
concerning the analysis of these field distributions are in- 
troduced in Sec lIVj more information on /iSR techniques 
in general can be found in Ref. [13. 



III. MAGNETIZATION STUDIES OF 
Lai.83Sro.irCu04 

The magnetization measurements in an applied mag- 
netic field of 20 mT are presented in Fig. [TJ The 
overall diamagnetic signals exhibit certain peculiarities. 
While the slowly cooled SFCW data represent the mag- 
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FIG. 1: (color online) Field-cooling magnetization curves of 
Lai.83Sro.i7Cu04 for ^loH — 20 mT and H || c for the differ- 
ent cooling procedures described in the text. 



netization curve of a rather equilibrated mixed state, 
the fast cooled FFCW measurement shows an initial 
low-temperature magnetization which deviates from the 
SFCW curve, reflecting a fairly undefined vortex con- 
figuration obtained by this cooling procedure. In con- 
trast, the SFCC data show a hump below T « 35 K. 
The field-cooled magnetization of type-II superconduc- 
tors during cooling (FCC) and during warming (FCW) 
were examined theoretically by Clem and Hao.^-'^ They 
demonstrated that the important parameters in analyz- 
ing the irreversible behavior of the low-field dc suscep- 
tibilities are the critical current density, the sample di- 
mensions, and the lower critical field Hd- However, the 
theoretical model developed in Ref. [Ill is subject to vari- 
ous simplifications and cannot be directly adapted to the 
curves presented in Fig. [T] Nevertheless, the detected 
hysteresis between the cooling and warming measure- 
ments is readily explained by differences in the flux expul- 
sion and reentering in the presence of vortex pinning.— 
In order to elucidate the origin of the distinct observa- 
tions in the case that the sample is cooled very fast, the 
magnetization relaxation of the FFCW and SFCC mixed 
states was studied over a time of about six hours for 
various temperatures between T = 5 K and as de- 
picted in Fig. [21 It is seen that for temperatures below 
about 10 K the vortex configuration hardly relaxes at all. 
For higher temperatures the magnetization returns to the 
"equilibrium one" with a maximum in the relaxation rate 
between about 20 K and 25 K which decreases again for 
temperatures above. Figure [3] shows the relaxation rate, 
defined as the slope of the measured magnetic moment 
in Fig. [2] for times t > 8100 s. The vortex relaxation 
for FFCW and SFCC is simflar for both studied situa- 
tions. The small observed differences in the relaxation 
rates seem to be merely related to the differing initial 
conditions. Even though a fully quantitative analysis of 
the vortex relaxation is beyond the scope of this article, 
the low-field flux-line dynamics appears to be different in 
at least the three temperature regions described before. 
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FIG. 2: (color online) Temporal evolution of the magnetic 
moment in an applied field of 20 mT in the mixed state 
of Lai.83Sro.i7Cu04 obtained for the different cooling pro- 
cedures, mo denotes the magnetic moment at f = defined 
for each temperature by data similar to those shown in Fig. [1] 



A similar behavior is found in Bi2Sr2CaCu208+5 where 
especially in a lovif magnetic field different vortex-pinning 
regimes have been identified. Yet, Lai.83Sro.i7Cu04 
is far less anisotropic than Bi2Sr2CaCu208+(5 and the 
more three-dimensional vortex structure is dominated 
by Josephson coupling rather than electromagnetic cou- 
pling. Thus, the pinning of individual "vortex pan- 
cakes" like in Bi2Sr2CaCu208+5 at low temperatures ap- 
pears unlikely. However, the overall observed magnetiza- 
tion relaxation still suggests a crossover between distinct 
vortex-pinning scenarios also in Lai.83Sro.i7Cu04. 
ft should be noted as well that the applied field of only 

lie 

20 mT at low temperatures is smaller than /io^ci in 
almost optimally doped Lai.83Sro.i7Cu04 As can be 
seen in Fig. [1] this does not affect the macroscopic "equi- 
librium magnetization", yet, the complete loss of dy- 
namics for T < 10 K might reflect the "final freeze-in" 
of the flux distribution at about the temperature where 
fioH^}^ w 20 mT." 

Similar differences between FCC and FCW magneti- 
zation data have also been observed for other applied 
fields (with slightly shifted temperature intervals). How- 
ever, previous /iSR experiments on a single crystal of 
Lai.83Sro.i7CuO^ have shown that the field of 20 mT 
is high enough to introduce a sufficient amount of vor- 
tices in the material which generate the characteristic 
magnetic-field distribution in the type-H superconduc- 
tor. Therefore, the magnitude of the applied field has 
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FIG. 3: (color online) Magnetic-moment relaxation rates in 
Lai.83Sro.i7Cu04 in fioH — 20 mT as a function of temper- 
ature for the two studied cooling procedures. The rates are 
determined from the data shown in Fig. [2] for t > 8100 s. 



also been chosen to be 20 mT for the present study. In 
the following it shall be investigated how the /zSR re- 
sults change for the various vortex configurations studied 
above. 



IV. ^SR STUDIES OF Lai gaSro 17CUO4 

For the initial /iSR studies the sample was cooled 
from above Tc to the lowest temperature in the applied 
field H II c at a cooling rate of w —25 K/min and the 
data were recorded while warming up the sample step- 
wise to T > Tc. This temperature cycle corresponds 
to the FFCW procedure of the magnetization measure- 
ments. First, the resulting data were analyzed using a 
sum of three oscillating signals with Gaussian damping 
( "3-Gaussian" method) which has been shown to be a 
reasonable approximation of the data if the probed field 
distribution originates from an ordered (or weakly dis- 
torted) vortex lattice.-^ In this case, the muon-decay 
asymmetry which is proportional to the spin polarization 
is modeled by 
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Here the Ai are the partial asymmetries of the three sig- 
nal contributions, ai the respective Gaussian depolariza- 
tion rates, and Bi the average fields of each of the frac- 
tions. The Ai are common parameters for all tempera- 
tures whereas the Ui and Bi are temperature-dependent. 
(f is the common initial phase of the muon spins with 
respect to the positron detector. The total decay asym- 
metry is ^ = Ai^ A2+ Ay, and the second central moment 
of the corresponding field distribution is given byi^ 



(2) 
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where 



1 ^ 



(3) 



is the first moment of the field distribution. Assum- 
ing the depolarization of the muon-spin ensemble is 
caused only by the inhomogeneous field distribution gen- 
erated by the vortex lattice and random nuclear mo- 
ments, then the contribution of the flux-line lattice to 
the second central moment is given by a^^ — a"^ — a"^, 
where (Tq = 0.23(1) /^s^^ is the muon-spin depolarization 
rate above Tc ~ 37 K in the normal state of the mate- 
rial. In order to further characterize the obtained field 
distributions also the dimensionless skewness parameter 

a = (AB^y^^ / (AB'^y^'^ is calculated; it represents 
the lineshape asymmetry of P(i3)i^ii£ Given the sum of 
Gaussian distributions, the third central moment can be 
written as 



{AB') = ^J2^^ - (S» 3 <J^hl + [B, {B)) 



A 



(4) 



For the final interpretations within the above model, the 
assumptioE^ii CTsc oc X^^ {\ab is the in-plane magnetic 
penetration depth) is essential. In order to test this as- 
sumption the data were also analyzed using an analytic 
Ginzburg-Landau (AGL) model J^ii^ Here the spatial dis- 
tribution of the magnetic field in the vortex state of a 
superconductor is modeled by the Fourier series 



B{v) = {B)Y,' 



— if 
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exp (— jKr) 



(5) 

where r = {x,y), K are the reciprocal lattice vectors of 
a two-dimensional hexagonal vortex lattice, K = |K|, 
and Ki is a modified Bessel function of the second kind. 
The parameter f^a representing the suppression of the 
superconducting order parameter due to overlapping vor- 
tex cores has been set equal to 1 (no suppression) in the 
analysis since the vortices are far apart in the low applied 
field used. The length is the effective vortex-core ra- 
dius which is a variable parameter within this model. The 
field distribution P{B) probed by the muons is obtained 
by random sampling of B(r) and finally the muon-decay 
asymmetry is given by 

A{t) = A exp {-(Tlt^/2) J P{B) cos (7^Si + Lp) AB. 

(6) 

The Gaussian prefactor in Eq. ^ accounts for broad- 
ening of the field distribution by nuclear dipole fields as 
well as weak random pinning^^ For the present data tig 
is growing continuously from 0.45 /is~^ just below Tc to 
0.55 Ais-i at T = 1.6 K (cf. Fig.gjl). 
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FIG. 4: (color online) Results of the analysis of the TP /^SR 
data on Lai.83Sro.i7Cu04: Temperature dependence of (a) 
obtained from the AGL model and g^c from the "3- 
Gaussian" analysis; the solid red line is a fit to the semi- 
classical model of Ref. [3^ using an order parameter with 
d^2_y2 symmetry, (b) the probed mean field, (c) the skew- 
ness parameter a, (d) the width of the Gaussian in Eq. ^ 
accounting for symmetric broadening of P{B) within the AGL 
model and the effective vortex-core radius ^v, and (e) the in- 
verse width of the Lorentzian weighting function in Eq. ^ 
given in units of the inverse inter-vortex distance of a fully 

ordered hexagonal flux- line lattice, a — "I>o/ (\/3 (B)) , 

where $o = h/{2e) = 2.07 x 10"^^ Tm^ is the flux quantum. 
For further details, see text. 
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The comparison of asc obtained by the "3-Gaussian" 
analysis (black squares in Fig. with A"^^ as deter- 
mined by the AGL analysis (blue diamonds in Fig. |4^) 
reveals that for this set of data both parameters cannot 
be simply scaled to match in the full temperature range, 
yet, both curves show pronounced changes in their 
curvatures around T w 13 K. The skewness parameter 
a is depicted in Fig. HJ:. Below T « 20 K a is virtually 
constant; the small differences in the absolute value 
of a between the various models are merely the result 
of effectively cutting off high-field contributions in the 
"3-Gaussian" analysis, thus leading to an overall smaller 
a which is less susceptible to minor effects in the "high- 
field tail" of P{B). One may note that the "3-Gaussian" 
analysis exhibits a strong drop and a sign change in a at 
T = 33(2) K. While such behavior usually is attributed 
to a melting of the vortex lattice,— here the detailed 
temperature dependence of a is most likely an artifact 
of the fitting procedure involving three Gaussians for an 
almost symmetric field distribution slightly below Tc- 
This view is supported by the fact that the irreversibility 
line in La2-a:Sr2;Cu04 at this low applied magnetic field 
tends to be much closer to Tcj ^^'^^ 

However, it is most important to check the temperature 
dependence of {B) determined by the different analyses. 
Figure |3}d shows {B) for the above described models 
as well as its values obtained by a maximum-entropy 
approacb^^ which has the advantage that it is not tied 
to any model. For the maximum-entropy analysis an 
apodization with a Gaussian (cTapod = 1 fJ-s) was applied 
to the asymmetry spectra. This leads to additional 
symmetric broadening of the determined field distribu- 
tion but does not change {B). The various analyses all 
yield qualitatively the same results: while the mean field 
is slightly diamagnetically shifted and approximately 
constant for T > 13 K, for T 13 K (B) strongly 
rises with decreasing temperature — even substantially 
above the applied field. Additionally, as shown in 
Fig. UJi, in order to explain the data obtained from 
the AGL fit would have to increase drastically in the 
same temperature range (cf. Fig. HJi) whereas it assumes 
values between about 10 nm and 15 nm above T « 13 K 
before it diverges close to Tc. 

Since the magnetization measurements have revealed 
a change in the vortex dynamics at T w 10 K and 
especially that the mixed state, generated in this initial 
/iSR experiment, is far from being equilibrated it shall 
further be investigated if the increase in (_B), ^v, and 
the second central moment of the fiSK field distributions 
at low temperatures is related to those observations. 
For this purpose, additional fiSK measurements have 
been performed. In a first step the single crystal has 
been cooled slowly in the applied magnetic field from 
above Tc to finally T = 1.6 K with a cooling rate of 
—0.2 K/min while the measurements have been done 
at intermediate stable temperature steps lasting about 
one hour each (SFCC). Subsequently, the so generated 
"equilibrium vortex state" served as starting point 
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FIG. 5: (color online) Normalized magnetic field distributions 
for selected temperatures obtained by a maximum-entropy 
analysis of the TF /iSR data on Lai,83Sro.i7Cu04 for the dif- 
ferent cooling procedures. The insets show the correspond- 
ing asymmetry spectra for T = 47 K (black squares) and 
T = 1.6 K (orange circles), respectively. The soUd lines in 
the insets are fits to the data using Eq. ((6)1. For T = 1.6 K, 
P{B) is calculated using Eqs. © and d?]); for T = 47 K > Tc, 
P{B) = S{B^fioH). 



for another series of measurements which have been 
conducted during the stepwise warming to above the 
critical temperature (SFCW). Figure [5] shows a com- 
parison between the static field distributions obtained 
by the maximum-entropy analysis of the TF /xSR data 
as well as selected asymmetry spectra for the different 
cooling procedures — the differences are only marginal. 
The further analysis with the above described models 
confirms this finding (cf. Figs. [S^-d and[7^-d). Also, 
the determined field distributions presented in Fig. [5] 
agree well with the expectations where, e. g. given the 
model parameters obtained by the AGL analysis at the 
lowest temperature the "high- field tail" of P{B) should 
range up to 32 mT. The strikingly similar results for 
the different vortex configurations indicate that the 
/iSR measurements at this small applied field mostly 
probe the local arrangements of the fiux lines, rather 
than their long-range order. This is especially seen 
in the broadening parameter dg which originally had 
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FIG. 6: (color online) Temperature dependence of the various 
model parameters for the SFCC TF fiSK measurements of 
Lai.83Sro.i7Cu04 — analogous to Fig. |4l 
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FIG. 7: (color online) Temperature dependence of the various 
model parameters for the SFCW TF /iSR measurements of 
Lai.83Sro.i7Cu04 — analogous to Fig. |4l 



been introduced in the model to take into account 
nuclear dipolar broadening and disorder of the vortex 
lattice. The temperature dependence of this parameter 
is almost identical for the different cooling procedures, 
even though the "vortex lattice" should be much more 
disordered and essentially not clearly defined in the 
initial FFCW measurement, thus strongly suggesting 
that the broadening of the field distribution is not just 
related to vortex-lattice disorder. 

The magnetization measurements do not yield any 
confirmation for an increased (B) due to the vortex 



lattice itself on the macroscopic scale — the applied field 
and the magnetization are constant below T 10 K 
(cf. Fig. [T|). Hence, it is unlikely that the sudden 
increase m (AB^) below T = 13 K— which is inevitably 
connected to the rising {B) — is related to an increase 
in superfluid density. Moreover, Aab(O) as obtained by 
the AGL model with the freely varying parameters is 
as low as 160 nm which seems to be far too short for 
Lai.83Sro.i7Cu04!^'^ Furthermore, while the observa- 
tion of the rather large effective vortex-core radius at 
13 K < T < Tc is in agreement with previous studies of 
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La2-a;Sr2;Cu04, where an expansion of the vortex cores 
has been found in a relatively low applied fieldjSS. the 
low-temperature values appear to be unreasonably large. 



V. POSSIBLE EXPLANATION FOR THE 
OBTAINED FIELD DISTRIBUTIONS 

The origin of the increase in the probed mean field, 
the second central moment of P{B), and the deduced 
effective vortex-core radius still remains to be clari- 
fied. As shown in Fig. |8] for the presently studied 
Lai.83Sro.i7Cu04 single crystal zero-field /xSR measure- 
ments indicate a slightly enhanced spin depolarization 
at low temperatures — in accordance with the fact that 
the so-called "cluster-spin-glass phase" in La2-2:Srj:Cu04 
might exist up to a doping level of x = 0.19.^^ Therefore, 
a small increase in the width of the probed field distribu- 
tion at low temperatures (T < 3 K) could be expected, 
however, static magnetic phases accounting for the dra- 
matic rise in {B) for T < 13 K can be excluded from 
these data. 

The temperature dependence of {B) is strongly reminis- 
cent of observations by /iSR in Bi2Sr2CaCu2 08-i-5j^'^ 
there in applied fields of 0.3 T to 0.4 T at low temper- 
atures the mean field probed by the muons also grows 
substantially with decreasing temperature whereas in 
the higher field of 1.5 T the probed {B) equals the 
applied field. It has been argued that in low fields 
when the vortex motion "freezes" at low temperatures, 
the random sampling of the spatial field distribution 
of the vortex lattice is perturbed since pinning cen- 
ters trap vortices and also offer possible muon stop- 
ping sites while at higher fields intervortex interactions 
dominate.— Of course, as already mentioned in Sec. IIIIl 
it has to be considered that the vortex-lattice topol- 
ogy is different in La2-xSra;Cu04 and the much more 
anisotropic Bi2Sr2CaCu2 08-i-5 system. The observed 
changes in the mean field in Bi2Sr2CaCu208+5 fall to- 
gether with the transition from a "vortex liquid" to a 
"vortex solid"— whereas the melting of the vortex lattice 
in the present sample is much closer to as discussed 
above already. Though, especially in relatively low mag- 
netic fields La2-xSr2:Cu04 shows a rich vortex-matter 
phase diagram where also in the solid phases the vortices 
are subject to thermal fiuctuations.'^*' Therefore, at high 
temperatures the vortices might "hop" thermally acti- 
vated between different pinning centers which causes the 
muons only to be exposed to an average field [thus leading 
to a broadened but still asymmetric P{B)] while at low 
temperatures these thermal fluctuations are reduced and 
more muon spins precess in the high static fields close to 
the vortex cores which have a substantial effective size in 
the very low applied magnetic field. This could explain 
the observed positive shifts in the first moment and the 
second central moment of the measured field distribu- 
tion. This view is supported by the magnetization data; 
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FIG. 8: (color online) Zero-field /iSR asymmetry spectra ob- 
tained for Lai.g3Sro.i7Cu04 at T = 47.4 K and T = 2.5 K. 
The low-temperature measurement shows a more exponential- 
like decay at small times and a slightly increased depolariza- 
tion rate. 



even though the relevant time scales differ in the ^SR 
and magnetization measurements, the observed almost 
cooling-cycle-independent change in the vortex dynam- 
ics around T — 10 K in the magnetization data as shown 
in Figs. [5] and [3] seems to be correlated with the pecu- 
liar variations in the field distributions probed by the 
muons. 

Very similar to the data discussed here — although in 
higher applied fields — Harshman et al. observed an in- 
fiection point in the temperature dependence of the sec- 
ond central moment of the field distribution measured 
by TF ^SR in a single crystal of YBa2Cu307_5.^ Also 
there, the increase of (AB^) at low T was attributed to 
vortex-pinning effects which hence have been taken into 
account in the analysis. This analysis, however, does not 
provide an explanation for the increase in {B) as shown 
above. Therefore, instead of introducing pinning effects 
into the model for Usc, in order to accomodate the pos- 
sible enhanced correlations between the muons and the 
vortex cores in the analysis, the random sampling of the 
magnetic fields within the flux-line lattice is replaced by 
a phenomenological weighting function w-e enhancing the 
contributions of the flelds close to the vortex cores to the 
resulting P{B). The two-dimensional spatial held distri- 
bution calculated using Eq. ([5]) is weighted by the (not 
normalized) Lorentzian 



1 



(7) 



where rj represents the distance to the center of the vor- 
tex core j and S is the inverse width of the Lorentzian. 
For the calculation three to six neighboring vortex cores 
are taken into account. The obtained widths of the 
weighting function are depicted in Figs.|4^,[6^, and[7^, re- 
spectively; above T = 20 K where P{B) is quite narrow, 
S is strongly correlated with dg but assumes generally 
small values, meaning that the flelds are almost sampled 
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randomly (2 = 0). For T ^ 20 K, 5 increases only mod- 
estly with decreasing temperature until T = 13 K be- 
low which temperature the vortex-core contributions to 
the overall field distribution grow intensely — consistent 
with the changes in P{B) observed during the previous 
analyses. On the other hand, {B) within this extended 
model remains essentially temperature-independent and 
slightly diamagnetically shifted with respect to the ap- 
plied field (cf. Figs. HId, IBh, and[7lD, red circles), overall 
consistent with the magnetization measurements. Due 
to the only small diamagnetic shift of {B) and the scat- 
tering of the determined mean-field values it cannot be 
finally decided if the observed difference in the absolute 
magnetic moment (cf. Fig. [1]) for the employed cooling 
cycles is reflected in {B). (Please note that {B) in this 
case only refers to the mean field connected to the pure 
randomly sampled ordered vortex lattice, whereas the 
overall {B) of course is not independent of temperature 
as discussed before.) Also, as seen in Figs.|3Ji,[BJi, andjTJi 
the deduced in this analysis is approximately constant 
as a function of temperature. The resulting tempera- 
ture dependence of A"^^ is shown in Figs. and 
[7^; it does not at all show pronounced changes of the 
curvature around T = 13 K and is approximately linear 
in T at low temperatures. The full temperature depen- 
dence of A~jf can be fairly described by the semi-classical 
model of Chandrasekhar and Einzelj^i^ where the pre- 
sented data are consistent with a weak-coupling dx2-y2 
order parameter with a maximum zero-temperature gap 
value 2 Ad(0) = 15.8(3) meV = 4.96(9) fceTc and a zero- 
temperature penetration depth Xab{0) = 199(3) nm (solid 
red line in Figs. H^, [6^, and [7^). The cited uncertain- 
ties are of statistical nature only. Especially the value 
of Aob(O) is subject to sizeable relative systematic errors 
of at least 5 % due to the phenomenological treatment 
of the vortex-pinning effects during the analysis. Nev- 
ertheless, the so determined value of Aa6(0) w 200 nm 
agrees well with data found in the literature^"''^^ while 
it appears to be estimated much too short if a randomly 
sampled vortex lattice was assumed as discussed at the 
end of Sec. |IVl 

It should be noted that the presented explanation for the 
observed field distributions might not be the only pos- 
sible one — yet, the introduced phenomenological model 
describes all obtained data consistently even though up 
to now the weight S cannot be fully related to the under- 
lying quantities such as the pinning potential. Any other 
model would also have to be coherent and particularly it 
would have to provide an explanation for the tempera- 
ture dependence of (B). 

One such other possibility to explain the /iSR data would 
be the presence of field- or vortex-induced antiferromag- 
netic order which would not be visible in the macroscopic 
magnetization measurements. This phenomenon has 
been reported to be present in cuprate high-temperature 
superconductors in applied fields of a few tesla by various 
techniques f2i"J£ Yet, based on neutron-scattering results 
which indicate no ordered magnetic phase in low fields 



in La2-xSra;Cu04 with such high doping level as used in 
the present study, Sonier et al. concluded from their /iSR 
studies that only disordered field-induced magnetism is 
likely to occur in low applied fieldsi^ Therefore, the pres- 
ence of an ordered magnetic state can be excluded for the 
presently discussed data. 



VI. DISCUSSION AND CONCLUSIONS 

The TF /iSR data presented in this work show a 
striking resemblance with data obtained earlier by 
Khasanov et al. in fiSK experiments on single crystals 
of cuprate high-temperature superconductors. Their 
data have been interpreted as evidence for two distinct 
energy gaps in the quasiparticle excitation spectrum, 
e. g. of La1.83Sro.17 CUO4," YBaaCusOr-a,^^ and 
YBa2Cu408'^ Note that this is basically different 
from reports of slight deviations from a dx^-y'^ gap 
symmetry as measured, e. g. by phase-sensitive methods 
in YBa2Cu307_54S The conclusions of Ref. [HI are 
mainly based on the observation of an inflection point in 
the temperature dependence of the second central mo- 
ment of the magnetic field distribution P{B) obtained 
by TF /iSR measurements in the vortex state (like 
seen in Fig. H^), while the corresponding temperature 
dependence of the measured mean field has not been 
discussed at all. Our TF /xSR data on Lai.83Sro.i7Cu04 
are very similar, yet, the above presented analysis within 
the framework of conventional vortex-lattice-generated 
field distributions extended by vortex-pinning effects 
appears to describe the overall /iSR and magnetization 
data in a coherent way. This is also in line with the 
suppression of the low-temperature increase of the 
second central moment of P{B) in higher applied fields^ 
where inter-vortex interactions compete with vortex 
pinning and minimize its effects on the measured field 
distributions. 

And, while it seems to be a reasonable approach to 
model two distinct gaps if an infiection point in A^^(r) 
is found (e.g. dj.2_y2 + sil or dr^2_y2 + idx^j^), the ob- 
servation of the strong increase in the probed mean field 
(B) in the TF /xSR measurements which is not observed 
in the magnetization measurements, renders it unlikely 
that the observed full second central moment oi P{B) is 
a true measure of A~^ at low temperatures. Moreover, 
it is pointed out, that the reported agreement between 
the obtained magnitudes of two-gap contributions from 
TF /iSR and neutron crystal-field spectroscopy^^iii is 
probably fortuitous since fundamentally different models 
were employed to obtain those values. While in the 
analysis of the TF /iSR data the contributions to the 
superfluid density were treated to be additive {after the 
integrations over the Fermi surface) r^i^ in the analysis 
of the neutron data the individual gaps are added 
together (before the integration over the Fermi surface) 
yielding effectively one pseudogap with an anisotropy 
Aa/Afcl 7^ 1 and shifted nodes. Taking into account 
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such an anisotropic superconducting gap in the modehng 
of the supcrfluid density would resuh only in tiny overall 
changes from the 6,^2 _y2 scenario since firstly, only the 
modulus of the gap is relevant and secondly, the inte- 
gration over the full Fermi surface^^ averages out small 
effects. Therefore, in general, a small gap anisotropy in 
the sense of Refs. or|4^ cannot be excluded, however, 
its determination from magnetic-penetration-depth data 
is rather difficult and indirect. 

In conclusion, using low-field TF /iSR in the mixed 
state of a Lai.83Sro.i7Cu04 single crystal a strong in- 
crease in the probed mean field and the second central 
moment of the measured field distribution is found at 
low temperatures. By combining /iSR and magnetiza- 
tion measurements it is shown that these effects seem 
to be primarily related to vortex-pinning effects chang- 
ing the sampling of the spatial field distribution by the 
muons and that the data can be consistently described by 
taking into account a single energy gap in the quasiparti- 
cle excitation spectrum with d3;2_j,2 symmetry. However, 
small variations of the superconducting order parameter 
cannot be excluded on the basis of the presented /iSR 
data. Given this assessment, the extrapolated magnetic 
penetration depth Aab(O) « 200 nm agrees well with re- 
ports of earlier experiments. 

Furthermore, this article commemorates that TF /xSR — 



in the way it has been employed for the present study — 
is first and foremost a very sensitive technique for the 
determination of local static or dynamic magnetic field 
distributions but the muon is neither a direct probe of 
the superfluid density of a superconductor in the mixed 
state nor of its order-parameter symmetry as also has 
been realized earlier.— Nevertheless, these parameters of 
a type-II superconductor can be deduced from TF /xSR 
data in the vortex state if a coherent description of the 
determined field distributions is available. 
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